Deformation and Depinning of Superconducting Vortices from Artificial Defects: A 

Ginzburg-Landau Study 
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Using Ginzburg-Landau theory, we have performed detailed studies of vortices in the presence of 
artificial defect arrays, for a thin film geometry. We show that when a vortex approaches the vicinity 
of a defect, an abrupt transition occurs in which the vortex core develops a "string" extending to 
the defect boundary, while simultaneously the supercurrents and associated magnetic fiux spread 
out and engulf the defect. Current induced depinning of vortices is shown to be dominated by the 
core string distortion in typical experimental situations. Experimental consequences of this unusual 
depinning behavior are discussed. 

PACS numbers: 74.20.-z, 74.60. Ge, 74.76.-w, 68.55.Ln, 61.46.+W 



I. INTRODUCTION 

One of the most important properties of superconduc- 
tors is their ability to carry currents without dissipation, 
allowing them to generate large magnetic fields. Many 
superconductors allow fields to penetrate in bundles of 
quantized magnetic flux, with associated whirlpools of 
current known as vortices. When these vortices are mo- 
bile, they spoil the perfect conductivity that make su- 
perconductors so useful. The quest to increase the max- 
imum dissipationless current Jc that a superconductor 
may carry has thus fueled intense study of vortex pin- 
ning. 

In recent years, pinning environments of artificially 
fabricated nanoscale defect arrays have been developed 
in hopes of better understanding and improving the pin- 
ning properties of superconductors. Some of the earlier 
contributions have involved macroscopic measurements 
(e.g. Jc, magnetization) performed on periodic "antidot" 
arrays^"'^. The antidot regions contain material which is 
rendered non-superconducting. Pinning behavior may be 
studied in these periodic systems using various imaging 
techniques^ . 

Much of the theoretical work on vortex pinning in de- 
fect arrays has employed numerical studies to focus^'^ 
on the behavior of large collections of vortices under 
the influence of a driving force (supercurrent). Molecu- 
lar dynamics^'^" approaches in particular usually employ 
simplified pinning potentials, in part to make possible 
simulations of large numbers of vortices, but also because 
information about pinning potentials at the microscopic 
level is simply unavailable. A few studies^^'^^ have fo- 
cused on energy scales for pinning varying numbers of 
flux quanta to the defects, as well as defect-vortex po- 
tentials as derived from the London equation^^'^^. How- 
ever, the latter approach does not allow for variation of 
the Cooper pair density, and in particular cannot cor- 
rectly treat the vortex core. In the London approach, 
vortex cores are usually assumed to be rigid in shape, 
and interactions of vortices with their environments are 



determined by the core position as well as the distribu- 
tion of currents"'^^'^^. Our results demonstrate that the 
vortex core in fact deforms dramatically near an artifi- 
cial defect: when the vortex center is sufRciently close to 
the defect, a string of suppressed order parameter devel- 
ops from the vortex position to the pinning center edge, 
while the currents and magnetic flux spread out over a 
large area (see Fig. 2 below). 

Usually, experimental studies of depinning are con- 
ducted in the presence of a net current across the su- 
perconducting film. The current exerts a Lorentz force 
on the flux quantum, dislodging the vortex from the de- 
fect for currents exceeding the depinning threshold. We 
find that, in addition to this, the external current distorts 
the conservative part of the pinning force for a sufRciently 
large net current, leading to important changes in the or- 
der parameter, supercurrents, and depinning threshold. 
Specifically, we find for many typical experimental situ- 
ations that in the presence of a driving current, the core 
strings associated with specific dots will stretch to reach 
neighboring dots. At the depinning threshold, the cur- 
rent associated with a vortex unwraps from one pinning 
site and encircles its neighbor, without the order param- 
eter ever forming a compact core region. Thus, for many 
experimental situations, interstitial vortices never form 
when vortices are depinned, as is commonly assumed. 
The depinning process should instead be understood as 
one in which the vortex cores are elongated by the driv- 
ing current, eventually interacting with more than one 
pinning site at a time and allowing the associated vortex 
currents to hop from site to site. 



II. METHODS 

Our calculations focus on two dimensional arrays of 
artificial pinning centers in the form of holes in a bulk 
superconductor. Our goal is to find the lowest energy 
state of the system for specified locations of a supercon- 
ducting vortex; from this we can construct a pinning po- 
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tential. The appropriate description of the superconduct- 
ing state is in terms of Ginzburg-Landau theory, which 
focuses on a complex superconducting order parameter 
tp{r), for which lipf is proportional to the local density 
of superconducting electrons. Unlike the London theory, 
Ginzburg-Landau theory is valid at scales as small as the 
coherence length, ^. Written in terms of dimensionless 
variables, the Ginzburg-Landau energy functional is 



Egl 



(v/z 
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In Eq. 1, ^ is the vector potential, and the magnetic field 
B{f) = V X A. K = is the Ginzburg parameter, the 
ratio of the magnetic penetration depth A and the co- 
herence length. Since the experiments on artificial defect 
arrays typically involve thin film geometries, we focus 
on them in our treatment. In considering the thin film 
limit, it is convenient to integrate by parts and replace 
the electromagnetic contribution, | V x with the quan- 
tity J ■ A. As the thin film limit is approached, one finds 
that ^, A, and J vary only very slowly in the z direction. 
The integration in the z direction is well-approximated 
by multiplication by d, the thickness of the film, and we 
may write 
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In Eq. 2, we use the definitions <j) = di/j, j = dJ, and 
(following standard convention^^) Xe/j = 

In this work, wc report on results obtained for Ae/ / = 
64 for defects larger than the coherence length. As is 
typically the case for thin films, the new length scale 
Ae// is considerably larger than A. In fact, this choice of 
Ae// is deep enough into the high Ae// limit that, apart 
from scale factors in Egl and 4>, the results vary little as 
Ae// is increased. 

To analyze the behavior of the vortex near the pinning 
center, we employ a mean field approach in which one 
minimizes Egl for a fixed vortex location to find A, 
and the current j. Our strategy for calculating ip and 
A self-consistently involves first holding (j) fixed at some 
initial guess, and minimizing Egl with respect to A and 
B. Next, we fix A and B and minimize with respect to 
(f). These steps are iterated until changes in the variables 
become negligible. In the case of an external current, mi- 
nor adjustments are needed which will be discussed later. 
We implement this self-consistent approach numerically 
by dividing the unit cell into a fine lattice of small unit 
cells. In this discrete scheme, 4>{f) is replaced by with 
ij specifying a grid point on a square lattice, while Aij 
and jij are defined on nearest-neighbor links between the 
grid points. Derivatives in Eq. 1 are replaced by the cor- 
responding finite differences. To model the defect array 



as accurately as possible, one desires a fine grid; we find 
that with a 128 x 128 mesh our results are well-converged 
with respect to the discretization. 

To see how one minimizes Egl under the constraint 
of a specified vortex location, it is useful to write the 
current j in terms of the order parameter cf) and the vector 
potential A. By minimizing Eql with respect to the 
vector potential and employing a Maxwell equation one 
has 
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0*(V||A-A)0 + (/.(-V|,/i-i')0*] (3) 



In Eq. 3 we have used for the order parameter (j) = |(/)| e*^, 
and the gradient V|| = (9x,9y,0). The familiar fluxoid 
quantization condition^^ arises from the requirement that 
the order parameter be single valued, i.e. § V9 ■ ds = 
2nny with riy an integer. Hence, in terms of j and A, 



2-Kny{ij) 



J/\(f>n ■ds+ (h A-ds 



(4) 



J/101') 



ds + $B 



In the second part of Eq. 4, $b is the total magnetic flux 
passing through the area of the contour, which we con- 
veniently choose to be the small unit cell associated with 
the grid point ■ij, while ny(ij) is the total number of "flux- 
oid quanta" contained in the contour of integration^^. It 
is through n„ that the vortex location(s) in the full unit 
cell of the system may be fixed; Uy = except at the 
grid points where we wish to place a vortex, for which 
riv = 1. As discussed in the Appendix, one can easily cal- 
culate $5 and A self-consistently in the thin film limit 
from the current j. Hence, armed with knowledge of 
and some specified realization of n„(ij), one can solve for 
j and A via Eq. 4. Using the expression for the current 
given in Eq. 3, one obtains V0; inserting V9 and A into 
Eq. 1 yields an expression depending only on |(^| and Ae/ / 
which we minimize with respect to \4>\ to obtain the order 
parameter modulus. 

Our calculational method is easily generalized to in- 
clude a supercTirrent jo flowing across the system, allow- 
ing the effects of a depinning current to be probed. This 
is accomplished by introducing a constant planar vector 
potential do- do causes a shift in the current given by 
AJ = -l^pao. 

In the above method, wc have been careful to explicitly 
take into account the thin film geometry typical of exper- 
iments on nanoscale defect arrays. However, we note that 
this approach may easily be modified to handle the bulk 
case, in which the vortices resemble filaments of magnetic 
flux instead of pancake structures. We have performed a 
few such calculations, and find that for antidot systems 
the results obtained are quite similar to the large Ae/ / re- 
sults we report here. This may be understood in terms of 
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the effective magnetic penetration depth for a thin film, 
Xeff = which is typically much larger than the bulk 
value A^^. This means that the energy stored in the mag- 
netic field generated by the supercurrents is quite small, 
so that the fact that the field varies as one moves out of 
the plane has little impact on the state of the system. 
The resulting energy functional is thus nearly identical 
to the bulk three dimensional case in the large k limit, 
with columnar antidots and vortices. 



III. RESULTS 

A. Zero Current 

At large vortex-defect separations, the core has the 
usual compact structure with supercurrents localized 
about it. Fig. 1 presents a perspective plot of \4>\, as well 
as a vector plot of the currents. The distances shown are 
in units of the coherence length, ^. Our choice of a unit 
cell with side spanning 20 coherence lengths and an an- 
tidot 5^ in size is typical of many of the nanoscale arrays 
studied experimentally. As the flux quantum nears the 
defect edge, it eventually reaches a critical distance dc 
(typically several ^, with precise value depending on the 
dot size and shape), where there is a sudden dramatic 
change. Fig. 2 illustrates the situation after the transi- 
tion: the vortex core has developed a string extending 
from the flux quantum position to the defect edge; si- 
multaneously the current now encircles the vortex-defect 
pair, and the magnetic flux created by these currents 
spreads over a larger area. The string is energetically 
favorable because it allows the formerly dense current of 
the vortex to spread out (engulfing the defect in the pro- 
cess), thereby reducing the kinetic energy of the state. 




FIG. 1. A perspective plot of the order parameter modu- 
lus, 101 (a) and current image (b) in a periodic antidot array 
just prior to the transition. Currents are localized about the 
vortex core, which has a compact structure. 
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(b) X/^ 

FIG. 2. Order parameter perspective plot (a) and currents 
(b) in the periodic array just after the transition. The currents 
circulate about the vortex-defect pair, and the vortex core has 
a string extending all the way to the defect edge. 



When a net supercurrent flows across a superconduc- 
tor, a Lorentz force is exerted on vortices; if these move 
in response energy is dissipated. In an artificial defect 
array, the driving force can be balanced by a pinning 
force given by the gradient of Egl with respect to the 



vortex position. Thus, one may regard Eql as the pin- 
ning potential. Fig. 3 shows this as a function of distance 
from the unit cell edge, which can be seen to have three 
distinct regions. For large separations, where the vortex 
core is compact, the pinning potential decreases relatively 
slowly. As the separation decreases it eventually crosses 
dc and the core-string structure appears. Hysteresis in 
the calculations, indicated with arrows, suggests that the 
transition is first order. However, as will be seen below, 
a finite current flowing across the unit cell changes the 
character of the transition, rendering it continuous. Be- 
low the transition one observes nearly linear behavior of 
the pinning potential, suggesting that the string carries 
an energy proportional to its length. The third region is 
announced by a discontinuous jump as the vortex is ab- 
sorbed by the defect, followed by a perfectly flat region 
inside the defect. This abrupt jump is due to a sud- 
den transformation of the order parameter. Immediately 
prior to absorption, constant Cooper pair density con- 
tours near the flux quantum have a semicircular profile. 
As the vortex enters the defect, this "semi-vortex" van- 
ishes, removing a finite amount of energy for an infinites- 
imal change in the vortex position. In later discussion, 
the Aej shall be used to refer to this jump in energy. 

B. Finite Current 

When a net current flows across the superconducting 
film, a Lorentz force is exerted on the current distribu- 
tion by the external magnetic field. This nonconservative 
force, given by -^o x J"'^^: acts in opposition to VEqj^. 
In terms of our units, the magnitude of the Lorentz force 
is Anj , where j is the net current magnitude (In CGS 

units, J = )■ Egl also depends on j . As will be 

seen, this dependence is particularly manifest near the 
depinning threshold. 

Figure 4 indicates the evolution of Eql as the current 
is increased. As the two graphs in Fig. 4 reveal, even 
for currents a reasonable fraction (17% and 33% respec- 
tively) of the bulk critical current Jc, Egl is similar to 
that of the zero-current case, except that the breadth of 
the linear region increases with drive current. Ultimately, 
as can be seen in the potential curve corresponding to the 
larger current, the linear behavior extends to either edge 
of the unit cell. This suggests there is a range of net cur- 
rent in which the vortex core always exists as a stringlike 
structure. Another important feature is the fact that the 
jump in energy when the fluxoid exits the dot (which we 
call the "ejection barrier" Aej) does not appreciably di- 
minish as the drive current j is increased. An interesting 
observation, apparent in Figs. 4, is that the Lorentz force 
increases with increasing current while the slope of the 
conservative part of the pinning potential remains nearly 
constant. Eventually, the Lorentz force exceeds the slope 
of Egl vs. fluxoid position before the jump Aej has been 
eliminated. This suggests that the ejection of the vortex 
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with increasing drive current will be sudden at low tem- 
peratures, i.e., once the flux quantum is able to exit the 
dot, the string tension will not be sufficiently strong to 
keep the vortex bound. We will see below that this is 
indeed the case. At higher temperatures, when T > A^j, 
one expects to see depinning for somewhat lower drive 
currents. 

Depinning occurs when the net current is made large 
enough to eliminate Agj. This occurs suddenly as one 
increases the current beyond a certain threshold, j^^. As 
Figure 5 indicates, the net supercurrent drops abruptly; 
presumably the remainder of the drive current appears 
in resistive channels. The fact that the solutions to the 
Ginzburg-Landau theory have a maximum possible value 
j even for large is highly analogous to what happens 
in thin superconducting wires and films^^. As in those 
cases, this maximum current should be identified as a 
critical current; in our case, it is the depinning current 
jdp- illustrated in Figure 6 and Figure 7, the depin- 
ning transition is marked by the sudden emergence of a 
core string. The string, a normal region which blocks 
supercurrents, accounts for the sudden drop in the su- 
percurrent flow across the unit cell. For the system pa- 
rameters we have chosen, the depinning current, is 

62% of the bulk critical current j^. 

Strikingly, currents capable of depinning the vortex do 
so without forming a flux quantum that has a compact 
core. In fact, as the vortex moves across the unit cell, its 
core always has a stringlike structure. As the flux quan- 
tum is driven from the pinning defect, a core "string" 
forms connecting the vortex to the edge of the pinning 
defect. Even before the flux quantum emerges from the 
pinning defect, this string is already partly formed. As 
the vortex traverses the unit cell toward the neighboring 
defect, the core string lengthens. This continues until the 
flux quantum reaches the unit cell edge, where the mag- 
netic flux and currents unwrap from the pinning defect 
and engulf the adjacent defect. Figure 8 and Figure 9 
illustrate this transition. Figure 8 displays the order pa- 
rameter modulus and current patterns just prior to the 
transition, while Figure 9 portrays the system just after 
the currents and flux have surrounded the neighboring 
defect. 

We conclude this section by noting that the depinning 
process found for typical antidot arrays is markedly dif- 
ferent than previous expectations, which commonly as- 
sume that in depinning the vortex takes a form simi- 
lar to Figure 1. We expect the latter scenario to hold 
when the interdot separation is very large compared to 
the coherence length ^. The different forms of depinning 
might be detectable in voltage noise experiments, which 
are analogous to shot noise experiments for electrons in 
conductors. When the vortices are small and compact, 
it is natural to suppose they will act as particles and, 
when depinned, yield a broadband noise proportional to 
the voltage across the system. If one increases the tem- 
perature and corresponding coherence length, the string 



depinning behavior described above must eventually ap- 
ply. This eliminates short wavelength and presumably 
high frequency dynamics from the system. We thus ex- 
pect a suppression of high frequency noise to accompany 
the transition into string-like depinning. 



Distance from unit cell edge (unit of £,) 

FIG. 3. Plot of the pinning potential showing broad linear 
regions and a gap in the potential at the defect boundary. 
The horizontal axis measures the distance of the vortex from 
the edge of the unit cell. Arrows indicate hysteresis in the 
calculation, a hallmark of a first order transition. 
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Distance from unit cell edge (a) constant component of vector potential, a° 



J =-.002 

X 




(b) magnitude of net current xio ' 



FIG. 4. Pinning potentials for = -.001 (a) and piG. 5. Plot of net current as a function of al (top) and 

Jx = -.002 (b) The exit barrier for the flux quantum does Ginzburg-Landau energy versus the current (bottom) 
not diminish despite increasing drive current. The slope of 
the dark line in both plots is the Lorentz force. 
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FIG. 6. A perspective plot of the order parameter modulus, pjQ 7 a perspective plot of the order parameter modulus, 

\(l>\ (a) and current image (b) in a periodic antidot array just |0| (a) and current image (b) in a periodic antidot array with 
prior to the transition. A core string has not appeared. core string, just beyond depinning 
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(b) X/^ 

FIG. 8. A perspective plot of the order parameter modulus, 
101 (a) and current image (b) in a periodic antidot array just 
prior to the transition. The core string is directed away from 
the pinning defect 



5 10 15 

(b) X/^ 

FIG. 9. A perspective plot of the order parameter modulus, 
101 (a) and current image (b) in a periodic antidot array just 
after the "unwrapping" transition. The core string is attached 
to the adjacent "destination" defect. Currents have moved 
from the pinning defect and encircle the "destination" defect. 
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IV. EXPERIMENTAL IMPLICATIONS 

Beyond the voltage noise behavior described above, the 
profile of the pinning potential shown in Fig. 3 suggests a 
set of measurements one might perform to seek an exper- 
imental signature of the unique pinning phenomena dis- 
cussed above. Past the depinning threshold, core strings 
moving among the antidots form linear channels which 
might appear as an anisotropy signal in neutron scat- 
tering measurements. Another possible test would in- 
volve working at temperatures large compared to Aej, 
but small compared to the energy at the unwrapping 
transition, as in Fig. 2. In this situation, an AC driv- 
ing force on the vortices, of magnitude small enough to 
leave them pinned by the strings (as might be provided in 
a vibrating reed experiments^) would allow the vortices 
to "rattle around" in the linear part of the pinning well 
illustrated in Fig. 3 and produce losses, whereas a DC 
force of equal magnitude would be dissipationless. Upon 
lowering of the temperature the vortices would be cap- 
tured by the defects, leading to lossless supercurrents for 
both AC and DC driving forces. An observation of these 
effects would yield indirect confirmation of the form of 
the pinning potential we find. Something like this may 
recently have been observed^''. 

Finally, a consequence of our results is that the de- 
pinning current should increase sharply as T is lowered 
below Aej. For T > Agj, vortices will be thermally ex- 
cited above the ejection barrier, and depinning can occur 
when the slope of Eql versus fluxoid position is equal 
to the Lorentz force. For T <^ A^j the current must be 
large enough to eliminate the ejection barrier to depin 
the vortex^^. 

Another interesting possibility is that a unique thermal 
depinning may occur as the temperature is increased in 
the regime of linear pinning. The presence of a string 
suggests that this may carry entropy at finite temper- 
ature, much as is the case of polymers. This entropy is 
proportional to the string length and temperature, and at 
high enough temperatures may overwhelm the energy per 
unit length found in our mean-field calculations. In anal- 
ogy with polymer behavior^*, this leads to unbounded 
growth of the string and effective depinning of the vor- 
tex. However, it is not clear whether the string remains 
sufficiently well-defined at the temperatures necessary for 
proliferation that the polymer analogy remains valid up 
to the transition. Further research into this possibility is 
currently underway. 

V. SUMMARY 

Using Ginzburg-Landau theory, we have given a de- 
tailed treatment of the microscopic aspects of pinning 
phenomena in nanoscale periodic arrays. Strikingly, we 
see an apparent first order transition involving the cre- 
ation of a string connecting the vortices to the defect, 



and an accompanying abrupt transformation of the su- 
percurrent and the magnetic field it generates. The string 
configuration leads to a region of linear pinning. Absorp- 
tion of the vortex by the antidot is marked by a jump in 
the pinning potential, Aej. Depinning occurs abruptly 
when the drive current exceeds Jdp- Various aspects of 
the pinning potential should be observable in experiment. 
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APPENDIX A: CALCULATION OF $s 



In this Appendix, we outline the calculation of the 
magnetic flux density in terms of the supercurrents in 
the thin film. The lattice variables mentioned in this 
work are discrete approximations of continuum quanti- 
ties. For the sake of brevity, however, the analysis in this 
Appendix is given in the continuum limit. The calcula- 
tion sketched here is largely parallel to an earlier work in 
the context of the London theory 

We adopt units for which the current J' satisfies J' = 
V X B. Taking advantage of the fact that V • A = in 
the London gauge, we find that 

J' = -V^l (Al) 

Fourier transforming Eq. Al yields 

= -fc^ik (A2) 

Defining fcjj = (fc^ + ^y)) have 

j'^:'' = -{kl + k^^)A^^y (A3) 

To express in terms of the two dimensional currents 

fiowing in the plane of the film one uses the fact that J' 
is essentially constant in the z direction in the thin film 
limit. Hence, we have 

jx,y^\ ,-2^Z^2' K (A4) 

1^ elsewhere J 
As the film thickness diminishes to zero, we obtain 

/-^U^o = <5(z)j// (A5) 
Fourier transforms equation A5 yields 

In Eq. A6 above, ky denotes the plane wave vector. We 
infer from Eqs. A3, A6 that 
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(A7) 



To obtain A^'^|^=o, and B^j^^o, an inverse Fourier trans- 
formation of Eq. A7 is needed. One finds that 



\z=0 



LxLy 



2{27rfJJ ' 



-(27r)-X 

fell 



-d^k (A8) 



(A9) 



To obtain B^\z=oj we use the fact that 

d 



2=0 



d_ 

dy 



\z=0 = 



2{2n) 



fell 



(AlO) 



In Fourier space, one has 

5k,| |.=o = -kyAl^, \z=o + fcxA^,, U=o (All) 
Hence, we find that 



dkxdky. (A12) 
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